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1. Introduction 

Accurate measurements of the muon properties were providing substantial information for 
the development of the elementary particle physics during many years. Nowadays precision 
experiments with muons serve as one of the basements of the Standard Model (SM) and 
give a possibility to look for new physics |fl, |2| . 

In this paper we discuss the present theoretical precision of the polarized muon decay 
spectrum description. The study is motivated by the experiment TWIST ||, which 
is currently running at Canada's National Laboratory TRIUMF. The experiment is going 
to measure the spectrum with the accuracy level of about 1 • 10 -4 . That will make a 
serious test of the space-time structure of the weak interaction. The experiment is able 
to put stringent limits on a bunch of parameters in models beyond SM, e.g., on the mass 
and the mixing angle of a possible right-handed VF-boson. To confront the experimental 
results with SM, adequately accurate theoretical predictions should be provided. This 
necessarily requires to calculate radiative corrections within the perturbative Quantum 
Electrodynamics (QED). In this paper we will concentrate on the effect of higher order 
leading and next-to-leading logarithmic corrections. 
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2. Preliminaries 



Within the Standard Model, the differential distribution of electrons (averaged over electron 
spin states) in the polarized muon decay reads 

_ j - j __r„(F W ±eP ( ,G(,)), r„ = ^, 

2mnE e 2rriam e 

C = COS0, X= — jj-^ 7J, Xq < X < 1, X = — o - ^ o) ( 2 -l) 



where and m e are the muon and electron masses; Gf is the Fermi coupling constant; 
9 is the angle between the muon polarization vector and the electron (or positron) 
momentum; E e and x are the energy and the energy fraction of e . Here we adopt the 
definition of the Fermi coupling constant following Ref. J?]]. Within the Standard Model, 
the muon decay happens due to the weak interaction of leptons and VF-bosons. The Fermi 
model corresponds to the limiting case of the infinite TV-boson mass. If Gf is defined 
according to Refs. ||, 0, the first order effect in the muon and VF-boson mass ratio gives 

3 ml \ 

1 + (2.2) 
5 m^J 

Functions F(x) and G(x) describe the isotropic and anisotropic parts of the spectrum, 
respectively. Within perturbative QED, they can be expanded in series in the fine structure 
constant a: 

F(x) = / Bom (x) + ^h(x) + (£)f2{x) + (£) 3 /s(x) + O (a 4 ) , (2.3) 

and in the same way for G(x). The Born-level functions /Born and <?Born are well known, 
including small terms suppressed by the m e /m^ mass ratio Q: 



/Born(z) = QX 1 + — § J 1 - -f x(l - x) + -p(A X 2 - 3x - xg) + 7/X (l 




5Born(x) = -2x^ 1 + -§ 1 - -=§ 1 - X + -i 4t - 3 - — X 



(2.4) 



where /?, 77, £, and 5 are the so-called Michel parameters || [9|, which in the Standard 
Model take the following values: p = 3/4, 77 = 0, £ = 1, and 5 = 3/4. By fitting the values of 
the parameters from the experimental data and comparing them with the SM predictions, 
the TWIST experiment is going to look for effects of non-standard interactions. 

The first order corrections fi(x) and gi(x) are also known with the exact account of the 
dependence on the electron mass |T^ , |j~3[ . Starting from O (a), radiative corrections to 
the electron spectrum contain so-called mass singularities in the form of the large logarithm 
L = ln(m^/mg) ~ 10.66. As demonstrated in Ref. (see also Table [l] below), the terms 
of the order O (otL) give the bulk of the first order correction. This is our reason to 
look first for higher order terms enhanced by the large logarithm. These enhanced terms 
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(excluding the ones related to the running of the QED coupling constant) cancel out in the 
expression for the total muon decay width at any order in a in accord with the Kinoshita- 
Lee-Nauenberg theorem [jD], |i~6|| . 

The second order corrections to the total muon decay width were calculated in Ref. jl7|] . 
At this order for the differential decay spectrum, we know only the leading logarithmic (LL) 
corrections [14| and the isotropic part in the next-to-leading logarithmic (NLL) approxi- 



mation 18]. The corresponding anisotropic part will be given below. The third order LL 



corrections will be presented here as well. 



3. The Fragmentation Function Approach 



First, I will describe briefly the application of the renormalization group method to the 
calculation of the leading and next -to-leading radiative corrections to the polarized muon 
decay spectrum. For a detailed foundation of the procedure and an extended discussion 
see Ref. 



The factorizations theorems, proved for QCD p9f| , can be easily translated to QED. 
In particular, by means of factorization, one can represent the differential spectrum of 
electrons ELS £1 convolution: 

d 2 r , 4 ^ r dz d 2 f 

[x 



J=e,7 x 

where d 2 Tj/(dzdc) is the differential distribution a muon decay with production of a mass- 
less electron (j = e) or a photon (j = 7) with energy fraction z {z = 2Ej/m fi , where Ej is 
the energy of the relevant particle) . To subtract collinear singularities from the differential 
distributions, we will use here the MS factorization scheme [^] with the scale fij. The 
fragmentation function Vj(x/z, fj,f,m e ) describes the conversion of a massless parton j into 
a massive physical electron. 

The spectrum of the massless parton can be expanded in a perturbative series: 

1 d2f ^. ,^ ., ^- m 



Tn dzdc 



(z,c,m„, f ) = If (z, c) + (*, c) + (^) Af (z,c) + (3 .2) 



where ct(/j,f) is the MS renormalized coupling constant, which will be converted further into 
the traditional QED on-shell coupling constant a ~ 1/137.036. The lowest order spectrum 
of massless electrons is defined by 

A£\z,c) = f (z)±cP tM g (z), f (z) = z 2 (3-2z), g (z) = z 2 {\ - 2z). (3.3) 

Here and in what follows, the sign before c should be chosen according to the charge of the 
decaying muon (plus for decay and vice versa). The O (a) correction to the massless 
electron spectrum reads 

A^(z,c) = fP(z)±cP,gi 1 \z), (3.4) 
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2z 2 (l - 2z)hx 



1-z 1 



4z + -z 6 In 



m. 



6 3 /' /ij 

+ 2z 2 (l - 2z) ( ln 2 (l - z)- 4Li 2 (1 - 2) - lnzln(l - z) - 2 In 2 z 



11 4 

T ~ 3z 



17 

T 



+ I — - — - : - -z 2 + yz 3 ) Inf.! - : i + 



1 



6z 2 - 6z A In 2 



J + 3z + 7 -z 2 + 3z 3 . 
b b 



(3.5) 



For the auxiliary photon spectrum with collinear singularities subtracted according to the 
MS prescription, we have 



4 0) (*, 



c) = 0, A^(z,c) = fW(z)±cP^(z), 
1 1 2 



3 3z 3 
2 \ 



-z 



+ H ( in ^ + 111(1 



11 



1 



+ \ ^- ^)lnz-- + — + —z - -zT - -z'\ 



3 3z J 



3 3z 12 



(3.6) 



(3.7) 



The isotropic parts of the first order corrections to the auxiliary massless parton distri- 
butions, fe(z) and f!y\z), are given by Eqs. (7,8) in Ref. |l8[| . By the choice of the 



factorization parameter value, fj,f ~ m», we avoid calculation of the unknown functions 

j (2) 

, since then they can not give rise to any LL or NLL correction. 

The fragmentation functions T>j can be obtained by solving the Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (DGLAP) evolution equations for QED, 



dVj(x, fj,f,m e 
dln/i 2 



J=e,f x 

where Pji are perturbative splitting functions, 

d(0) 



1 

22 / T ^«(M/)Pi(-,M/,m e J, (3.8) 



Pji{z,a(n f )) 



2tt 



(3.9) 



The initial conditions, which are required to solve the DGLAP equations by iterations, can 
be borrowed from QCD studies [21]: 



V™(x,v ,m e ) = 5(l-x) + ^-d 1 (x, f i ,m e ) + O(a 2 ) 



di(x,fj, ,m e 



2vr 

1 + ' V ln4 

mi 



1-x 



21n(l -x) - 1 



«Oo) , Mo 



ln^)(x) + 0(a 2 ). 



2ir mj 

The relevant lowest order splitting functions are 

'l + x 2 ] „ m/ , l + (l-x) 2 



(3.10) 
(3.11) 



1 - X 



^ 0) (x) 



p(°)(rc) = x 2 + (1-x) 2 . (3.12) 
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The plus prescription works as usually: 
l l 
dx [V(x)] + W(x) = J dx V(x)[W(x)Q(x - x min ) - W(l)), (3.13) 
o 

1 for x > 
for x < 



9(s) 



In a measurement of the muon decay spectrum, events with more than one electron in 
the final state require a special treatment. Starting from the second order in a, we have a 
certain contribution due to emission of real and virtual e + e~ pairs. Presumably a Monte 
Carlo event generator is needed to simulate the process of muon decay with pair production. 
Nevertheless, we will calculate the corresponding effect under a simple assumption, that 
an event with two electrons in the final state is treated as a pair of simultaneous muon 
decays. In order to have the possibility to drop the pair contributions (if they are taken 
into account in a Monte Carlo program), we decompose our results according to the classes 



of the corresponding Feynman diagrams in the same way as in Ref. [22|. Moreover, the 
decomposition will help us to demonstrate the cancellation of the mass singularities in the 
integrated decay width. Our results for pure pair corrections can serve further also as a 
benchmark for the Monte Carlo program. The next-to-leading electron splitting function 
can be decomposed into four parts: 



P£H*) = P&Xx) + ^ie 1,NS) (z) + J& 1,S) (z) + P£ M) W, (3-14) 



where Pil' 7 ^(x) is provided by the set of Feynman diagrams with pure photonic corrections; 
Pee' ^0^) is related to the corrections due to non— singlet real and virtual pairs; Pee 
stands for the singlet pair production contribution; and Pee' mt \x) describes the interfer- 
ence of real singlet and non-singlet pairs. By extracting the appropriate color structures 



from the known QCD results [23, 24, p^, 26 1, the explicit expressions for these functions 



were given in Ref. ||18[| . Here and in what follows, in the language of Feynman diagrams, 
the situation when the registered electron is connected by a solid fermion line with the gen- 
uine electroweak decay vertex is called non-singlet. The case, when the observed electron 
belongs to a pair produced via a virtual photon, is called singlet. 

The relation between the MS and the on-shell coupling constants reads 

, ^ a 2 . lA a 3 _ /4 15a 3 _ , 4N tn 

5 <"/>= a+ ^ ln 4 + 4? ln 4 + i^ + °( Q '>- (3 - 15) 

It is convenient to choose the renormalization scale to be 

Ho = m e . (3.16) 



Now we have everything for solving the DGLAP equations fl3.8|) . Using iterations for 
the electron fragmentation function decomposed into four parts, we get 

V e (x,Hf,m e ) = pW(x)+D( NS >(x)+pf(x)+pM(x), (3.17) 
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V™{x) = 8(1 -x) + ^d x (x,^m e ) + £-L f pW(a 



+ 7T- 



Q' 

27r y V2 



2vr iV c/ 2vr 



(3.18) 



a \ /I 



27r7 V3 

a > 3 
2vr 



^) 4(^ie 0) ^^ 0) (x) + |pi e °)(x)), 



(3.19) 



+ r- 



a 
2~^ 



Lj (^ 7 0) 55 PW 55 P<? (*) - ipW Pj? (x] 



L f = In 



(3.20) 
(3.21) 



where we systematically omitted terms of the following orders: O (a 2 L® \ , O (a^L 2 , 
O (ct 4 ), and higher. The photon fragmentation function at the lowest order, 



P 7 (x,/i /) m e ) = ^L f P^(x) + O (a 2 ) , 



(3.22) 



is sufficient for our purposes. The convolution operation is defined in the standard way: 

l l l 



dz ,„(x 
— A(z)B[ - 

z \z 



(3.23) 



A <g> B(x) = J dz J dz' 5(x - zz')A(z)B(z') = 

x 

The leading logarithmic terms for the QED fragmentation function are known up to the 
fifth order in a [p7f| . But, as will be seen from the numerical results, keeping contributions 
up to the third order is enough for the moment. 

The explicit expressions for the functions, which appear in the O (a 2 ) terms of the 
function T> e , are given in Ref. [18]. In the third order we have two more functions [28]: 



P e ( e 0) ® P e ( e 0) ® P e ( e 0) (x) = 5(1 - x) ( 16Cs - ^ ) ! 



24 



9 1 . 

H C2 

32 2 s 



+ 24 



1 + x ( 1 1 2 1 , , ,i ^ 
in x mxm(l — x) 



1-x V 12 



' ^Iln 2 (l-x) + ^ln!l - .,■) 



In a; 



l-i \2 



+6(1 +x)lnxln(l - x) - 12(1 - x) ln(l - x) + -(5 - 3x)lnx - 3(1 - x) 
3 

-~(1 + x) In 2 x + 6(1 + x)Li 2 (1 - x) , 
P®> ® Pj? ® P e ( e 0) (x) = (1 + x) ^4 ln(l - x) In x - In 2 x + 4Li 2 (1 - x) J 



2 /4 



+-(3x + 4x 2 )lnx + - - + 3-3x-4x 2 ln(l - x) (1-x) 



3 V x 



23 



(3.24) 



(3.25) 
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where 



oo 1 

T — 



k n 6 

k=l 



ti = 4, Li 2 (,).-/ dy Ml^). ( 3. 26 ) 



By convolution of the fragmentation functions [Eqs. ( |3.17| ) and ( 3.22 )] with the dif- 



ferential distributions ( |3,2| ), we receive higher order corrections to the electron spectrum, 
which will be presented in Sect. ||. In the results we can fix (see discussion in Ref. |H|| ) the 
factorization scale 

fj, f = m tJl (3.27) 

and get Lf — > L. 

4. Exponentiation 

Looking at the end point of the energy spectrum [x — » 1) of unpolarized muon decay, 
one can recognize that the first order correction becomes there negative and very large, 
making the result senseless. An extensive discussion of the phenomenon can be found in 
Refs. [12, 29]. The divergence is a clear signal to look beyond the first order approxima- 



tion. In fact, the Yennie-Frautschi-Suura theorem |30| allows to make a re-summation 
of the dangerous terms and to convert them into a definitely positive exponent. The ex- 
ponentiation procedure is not unique, it permits to involve also some terms convergent at 
x — > 1. In our case, the exponentiation is allowed to add several terms of the following 
types: O (a 2 L°) , O (a n L m ) with n > 3, < m < n, and O (a n L n ), n > 4. 

Let us consider two ways of exponentiation. In the first case one starts from the 
corrected cross section and tries to perform a re-summation of the known terms, which are 
divergent at x — > 1, by converting them into an exponent: 



F(x) 



x->l 



1 + — (L - 2) ln(l -x) + ... — * exp \—(L — 2) ln(l - x)\ . (4.1) 

7T I 7T J 



This is a kind of the so-called ad hoc exponentiation. The effect (see Table [l]) of this 
approach can be represented by the relative contribution of new terms generated by the 
exponent, 



6 e *l( x ) = exp{^(L-2)ln(l-x)} - 1 - -(£ - 2) ln(l 



x 



~ I Q V - 4L) ln 2 (l -x)-± {^L ln(l - x)\ \ (4.2) 

The most significant term above is of the order a 2 ln 2 (l — x) and gives a numerically 
important contribution at large x. Note that all the subtracted terms do appear in the 
perturbative results. 

The next step should be to check that in higher orders the exponent doesn't contradict 
the known (or anticipated) results. The above procedure in the case of muon decay can be 
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criticized [31], because the higher order leading logarithmic terms represent a mass singu- 
larity: one can not guarantee that all the large logarithms, coming from S^(x), disappear 
after the integration over the energy fraction. Nevertheless, the ad hoc exponentiation is 
not supposed to produce a complete result. The region of its application is limited: it deals 
with the terms, which are the most important in the end of the spectrum. 

There is another way of exponentiation, which avoids the problem of improper mass 
singularities. One can use the exponentiated representation of the electron structure (frag- 
mentation) function suggested in Ref. [32], which obeys the proper normalization: 

i 

j V™{x)6x = 1. (4.3) 
o 

The exponentiated structure functions are based on the exact solutions of the QED evolu- 
tion equations in the limiting case of soft radiation. For computations I used an advanced 
formula from Ref. p^l , where I substituted (L — 1), which was natural for e + e~ annihi- 
lation, by {L — 2). This substitution has not been obvious from the beginning, but it is 
clearly seen in the above ad hoc exponentiation. The usual {L — 1) factor corresponds to 
soft radiation off the electron, while the additional (—1) is due to radiation off the muon. 
In fact one can introduce the muon structure function into the master equation fl3.1|) , as 



discussed in Ref. [34|. The muon structure function does not give any large logarithms in 
our case. But still it can be used to describe the contribution of soft photon radiation off 
the muon. And the corresponding factor is just that (—1) instead of the usual (L — 1). We 
see, that it is easy to make a mistake in exponentiation of soft gluons (photons) in decay 
processes by forgetting about the radiation off the decaying particle. 

The relevant for us electron structure function |33| is taken within the leading logarithm 
approximation for pure photonic corrections with terms up to O (a 3 L 3 ) and supplied with 
exponentiation of some terms in higher orders. Convolution with the Born-level functions 
gives 

(*) = ® /o (*) , GgP (x) = Pg ® 90 (x) • (4.4) 

A subtraction of the known terms in the lower orders (n < 3) of the perturbative expansion, 
similar to Eq. (|4.2j), is used to receive the value of the relative contribution <5gp P (x). The 
latter contains also terms of the order O (a 2 L°Y which are not singular at x — > 1. This is 
due to the fact that soft radiation is allowed not only at the end of the energy spectrum, 
but in any other point in x. The resulting effect is spread all over the spectrum and appears 
to be visible at the two ends x — > and x — > 1. Numerical results (see Table ||) of the two 
approaches are close to each other in the large- x region (where an analytical agreement 
between the approaches can be observed as well). Contrary to the exponentiation of soft 



gluons in QCD [34], our procedure doesn't suffer from the renormalization scheme (and 
scale) dependence and from non-perturbative effects. 

Simultaneous exponentiation of photonic and pair corrections can be constructed as 
well. But it was criticized |35|, since soft pairs (contrary to soft photons) have a non-zero 
production threshold, which can't be taken into account by exponentiation properly. 
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5. Results 



To the best of our present knowledge, we can write now function F(x) from the master 
formula (^4|) as follows: 

2 



F(x) = f Bom (x) + —f 1 (x)+l — 



/ 2 ^(x )+ ^r s) (x) + /f s) (.) 



L 2 

2 



jf%) + /™(x) + ff(x) + /^(s) 
3r i L 3 



e(l,NS), 



c(l>S), 



,(l,int). 



, a 

+ f 27 



/^(x) + /f NS )(x) + /f S )(x) 



— + A/ Cxp (x) 



+ 0(a 2 L°,a 3 L 2 ,a 4 L 4 ). 



(5.1) 



Function G(x) takes the same form with the substitution /j — > ft. Effects due to virtual 
hadronic, and t + t~ pairs [36] are not shown explicitly, since they are of the order 

O (a 2 L°) . The effect of exponentiation is given by 



A/ exp (x) = f (x)5 s j(x), Ag cxp (x) = g (x)5 s j(x). 



(5.2) 



5.1 Analytical Results 

In the second order for the NLL corrections to the anisotropic part of the electron energy 
distribution, we have 

g}' i) {x) = 4x 2 (l - 2x) ^Li 3 (1 - x) + Si j2 (1 - z) - 2Li 2 (1 - x) m(l - x) 



+ lnxln (1 — x) — 3 In xln(l — x) + In x — £2 lnx + — £3 



14 



92 



+ I 6x + 24x x ) LI2 (1 — x) + I 6x — 5 x I In x ln(l — x 



86 



+ 
+ 
+ 



3x 

— - 12x + —x 2 + 8x 3 ) ln 2 (l - x) + ( — + 18x 2 - — x 3 ) In 2 x 
3x 3 / \ 12 3 



13 37 50 o 32 ., . , , 

1 x H x z x J ln(l - x) + 

3 3 3 3 1 



25 59 2 32 ., . . 

x + 6x H x lnx 

12 6 9 



8 29 9 o\ A 817 91 62 „ 607 , 

-8 H h 12x x 2 - 2x 3 C2 H a; H x 2 x 3 , 

3x 3 r 216 12 3 54 ' 



(5.3) 



^ 1,NS) (x) = 4x 2 (l - 2x) ( -Li 2 (1 - x) - ± lnxln(l - x) + \ ln 2 (l - x) 



1 



1 



/22 8 



- In 2 x C2 + 4x - 6x 2 + 12x d ln(l - x 



V 9 9x 



+ 



1 



76 



9 9 



9 



7 



86 



20 



— I — x x lnx 1 — xH x x, 



18 



9 



(5.4) 



\j\2 (1 — x) + In x ln(l — x) ) ( —x 2 



^ + \-x 



1 



+ 



1 

9 ~ 9x 



+ x + -x" 



9 y 



] ln(l - x) + 



3 2 
4 



- In x 



9 9x 



5 5 . . , 

H — x H — x mx 
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1 4 7 43 2 10 o 

H 1 x x z H x 6 , (5.5) 

3x 3 18 18 9 y J 

g£ M \x) = 4x 2 (l - 2x) ^Li 3 (1 - x) - 2Si, 2 (1 - x) - Li 2 (1 - x) lnxj 

+ (-\- Ux 2 + yx 3 j Li 2 (1 - x) + (-3x 2 + ^x 3 ) In 2 x 
1 1 28 2 \ 10 1 37 2 104 3 

— I — x x mxH x x H x . 5.6) 

333/ 9339 y 1 

The polylogarithm functions are denned as 

t- t \- f a Li 2 (y) a ^ 1 L ln2 ( 1 - y) ^ 

Li 3 (x) = / dy , Si, 2 {x) = - / dy . (5.7) 

J V 2 J y 

o o 

The O (a 2 L 2 ) corrections ff' j \x) and gf j \x) (j = 7, NS, S) can be found in Ref. @. 
Explicit expressions for the second order next-to-leading corrections to the isotropic part 
of the spectrum (f^ 1 (x), i = 7, NS, S, int) are given in Ref. [18|. 
The third order LL photonic contributions read 



/ 3 (0 ' 7) (x) = 8x 2 (3 - 2x)*(x) + (10 + 24x - 48x 2 + 32x 3 ) ln 2 (l - x 
_5_ 
12 



+ ( — + x - 8x 2 + 16X 3 ) ln z x + (-5 - 12x + 48x z - 64x a ) In x ln(l - x) 



+ (5 + 12x - 32x 3 )Li 2 (1 - x) + (-10 - 24x + 48x 2 - 32x 3 )C 2 

/ 13 21 64 3 \ , (11 _ „ 2 64 <A , m 
+ x H x 6 lnx + h 17x + 16x z x d ln(l - x) 

V 18 2 3 / V 6 3 / 

569 4 16 2 128 3 

H h -x x 2 H x 3 , (5.8) 

216 3 3 27 v ' 



(0>t)^ _ s Jn + (_ 2 _ 48x 2 + 32x 3 ) ln 2 ( 

8x 2 + 16x 3 ) In 2 x + (1 + 48x 2 - 64x d ) In x ln(l - x) 



gf J) (x) = 8x 2 (l - 2x)^(x) + (-2 - 48x 2 + 32x 3 ) ln 2 (l - x) 
1 

12 

+ (-1 - 32x 3 )Li 2 (1 - x) + (2 + 48x 2 - 32x 3 )C 2 

/ 5 5 64 .A / 7 o 64 „\ , , 

+ ( — + -x+ yx 3 j lnx+ ( ---7x + 16x 2 -yx 3 jln(l-x) 

133 13 16 , 128 o . . 

x x 2 H x 3 , (5.9) 

216 6 3 27 v 1 

1 3 

V(x) = 3Li 3 (1 - x) - 2Si 2 (1 - x) + ln 3 (l - x) - - In 3 x + - In 2 x ln(l - x) 

6 2 

- 3 In x ln 2 (l - x) - 3Li 2 (1 - x) ln(l - x) + 2C 3 - 3C 2 In — - . 

x 

And the third order LL pair corrections 1 are 

/(°. NS )(x) = 8x 2 (3 - 2x)$(x) + ( f + 16x - ^x 2 + ^x 3 ) ln(l - x) 



1 Strictly speaking, we have here pair and photonic corrections simultaneously. 
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+ 



= 8x 2 (: 



' 5 , 32 2 160 3 \ 73 67 16 , 128 „ , r 

-- - 4* + y x 2 - — x 3 j Inx + - + y x + y x 2 - —x 3 , (5.10) 

2/ v*/ x ( 4 272 2 160 ,\, , 

x 2 (l - 2x)$(x) + x 2 H x 3 ln(l - z) 

V 3 9 9 y 

' 1 ien \ 29 7 16 2 128 

x H x x' 

54 3 9 27 



1 128 2 160 ,\, 29 7 

— I x x mx . 

3 9 9 y 54 3 

" ' — 4x 2 In 2 x 



+ u 

= Q + 4x + 4x 2 ^) f 4Li 2 (: 

56 2 32 3 \ /n , / 29 i i 
H hl2x x 2 x 3 ln(l-x) + ™ 

3 9 y v ; V 9 3 

86 2 80 o 
— x z + — x d , 



(5-11) 



1 — x) + 41nxln(l — x) — In 2 x 

— x + 16x 2 

9 3 



/68 8 

+ U + 3^ 

32 ,\ , 287 4 

— x mx 

9 J 27 9x 9 

— x) + 4 In x ln(l — x) — In 2 x^ — — x 2 In 2 x 



+ 



! 13 cM) CM) 

-x H x H s 

9 27 



(5.12) 



8 2 32 o \ , . 
4x H — x x m(l — x) 

9 9 y V ' 



4 8 8 

h 4x + -a 

9 9x 9 

1 + 



+ V 9 9x 

16 2 32 o\ , 31 4 
- X + Y X J ln ^ + 27 + 27^ 



35 



oo iu 2 ou 

"Y X_ 27 X + 27' 



<3?(x) 



I 16 2 32 

+ 2x x 2 H 

\9 9 9 

= — In 2 x + ln 2 (l — x) — 21nxln(l — x) — Li 2 (1 — x) — C 2 . 



10 9 80 o . 
-x 2 + -x 3 ,(5.13) 



Functions <? 2 1,S ^ an d 92 > shown above, as well as / 2 1,S ^ and /.j 1 ' , which are given 
in Ref. IS], coincide with the results of my calculations starting directly from Feynman 
diagrams and using methods described in Ref. |}7]]. 



5.2 Cancellation of Mass Singularities 

An important check of the results is to demonstrate the cancellation of mass singularities 
in the total decay width. The cancellation of the mass singularities in the LL contributions 
due to photons and non-singlet pairs is rather trivial: 

1 1 

J dx f^\x) = j dx g^\x) =0, j = NS, 7, n = 1,2,3, . . . (5.14) 


It is guaranteed by the normalization conditions of the corresponding LL fragmentation 
functions. 

Now we should note that a naive integration of the electron spectrum gives rather the 
counting rate of electrons than the total muon decay width, since the number of the final 
state electrons can exceed the number of decaying muons because of real e + e~ pair emission. 
In other words, we should avoid the double counting of electrons in the contributions due 
to real pair emission. For this purpose, we can keep the non-singlet pair contributions and 
drop the singlet ones (see their definition on page ||). 
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Functions J^ 1 ' 111 *) and <?2 contain the double counting too. To resolve this problem 
we can use the splitting function [23, 38] 



1 + x 2 ( 1 
2 i ( -2Li 2 (-x) - 2 In x ln(l + x) + - In 2 x - C2 

2(1 + x)lnx + 4(1 - x), 



(5.15) 



which describes the transition of an electron into a positron in the relevant set of Feynman 
diagrams. The corresponding contribution can be constructed by convolution with the 
lowest order functions /o and go integrated over the positron energy fraction. One can see 
now that in the interference contribution the number of electrons is really twice as large as 
the number of positrons: 



|dx/ 2 (1 < int) (x) = 2 



/ dxP^ 



® fo(x) 



13 



3C 2 + 2Cs • 



(5.16) 



By resolving the problem of double counting in the contribution of the S-NS pair interfer- 
ence, we arrive to the cancellation of mass singularities in the following form: 



1 

/■ 



dx (/ 2 (1 ^(x) + ft im '{x) - P£™> ® /„(*)) 



f (l,int) 



(l,int) 



(5.17) 



Let us look now at the integral of the second order NLL non-singlet pair correction. It 
is known [|l7]] , that the integrated contribution of this correction contains large logarithms 
due to the running of the coupling constant: 



a 
2^ 



\jdxf^\x) 



Aa(m / , 
2T" 



dx fi(x) 



m e — >0 



(5.18) 



where 



Q- 



Aa(m u ) = — L = a(m„) - a + O (a 2 L u ) . 
3tt v 



(5.19) 



To demonstrate that using our results, note first that the relevant function consists of two 
parts: 



f^ S \x) = P^®f (x) + Z -f l (x). 



(5.20) 



One can check that 



dxP£ N ^ ®/ (x) = 0. 



(5.21) 
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It remains now to recognize that 



l l 



dx fx(x 
o " o 



dxMx), (5.22) 



which can be verified easily. Thus we checked successfully an important property of our 
analytical results. 

The anisotropic contributions to the decay spectrum can be treated in the same way. 
They don't contribute to the total decay width at all, because they vanish after the inte- 
gration over the angle. Nevertheless, the cancellation of mass singularities can be observed 
also in the forward-backward asymmetry of the decay, which is not affected by isotropic 
functions on the contrary. In particular, we have an equality analogous to Eq. (|5.16| ): 

l l 



dxg^'(x) = 2 J dxP^ 1 > ® 9o(x) = -- + &- (5.23) 



5.3 Numerical Results 

Now we can make numerical estimates of the effects due to higher order radiative correc- 
tions. We should confront them with the 1 • 10 -4 precision level of the TWTST experiment. 
In a typical experiment on muon decays one has almost 100% longitudinal polarization of 
muons (P„ = 1), since the muons are coming polarized from pion decays. In order to 
extract information about the Michel parameters £ and 5 one should study carefully the 
angular distribution of the electrons. But as can be seen from the analytical formulae, 
the dependence on the angle is rather simple and smooth. For this reason we restrict our 
numerical illustrations to the consideration of ^-dependence at several fixed values of the 
angle. 

Numerical results for the pure photonic corrections are presented in Table |], where 
we give the values of different contributions normalized by the Born distribution in the 
following way: 

s _ a_ fi(x) ±cP IM g 1 (x) 

2VT / B orn (x) ± cP^Born (x) ' 

"f { y\x)±cP,gt'\x) 





( a 


n\ 






a \ 


< 






2tt) 



fBom{x) ± cP^gBornix) 



n = 1,2,3 , 



4'-'=^-yf 7 ' ( 1 ±CP ^ 1,7 ' ( 1. (5.24) 

\2itJ / B orn(2;) ± cP M 5Born(^) 

One can see that the first order LL correction j| 0,7 ^ (look for f[° (x) and gf 1 (x) in 



Ref. [14]) provides the bulk of the effect, especially in the region of intermediate and large 
x-values. Convergence of the corresponding series in L in the second order corrections 
doesn't look so good: the NLL contribution is only about two times less than the LL one. 

There is a trick, which allows to make a better approximation in the region of small 
x. Looking closely at the argument of the large logarithm during the actual calculations of 
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X 


10 di 


10 o{ ' 


1 n4 r(0>7) 
10 t>2 


10 o 2 


10 O3 y 




1 r»4 c-exp 


0.05 


4590.1 


10325.0 


184.96 


-247.63 


-0.43 


0.00 


6.90 


0.1 


1715.1 


3257.7 


33.18 


-37.79 


-0.35 


0.00 


1.39 


0.2 


674.0 


1106.2 


-1.28 


0.34 


-0.15 


0.01 


0.05 


O.O 


364.0 


549.1 


—0.58 


/l AQ 

4.48 


—0.05 


0.01 


n on 
— 0.20 


0.5 


64.1 


82.6 


-6.50 


3.73 


0.05 


0.06 


-0.26 


0.7 


-160.3 


-214.9 


-1.28 


0.70 


0.07 


0.18 


-0.09 


0.9 


-470.3 


-592.1 


15.19 


-5.97 


-0.14 


0.72 


0.62 


0.99 


-971.9 


-1198.8 


69.84 


-26.08 


-1.95 


3.47 


3.70 


0.999 


-1439.8 


-1772.5 


155.10 


-57.86 


-6.64 


9.17 


9.75 



Table 1: Photonic corrections to [i decay spectrum versus x for c = 1 and P M = 1. 



integrals over the phase space of real photons, one can notice that it is rather x 2 m'f l /m i 



fit ""ei 

than simply m^/mj. So the modification L — » L + 21na; can be done in our formulae. This 
will move some terms from the sub-leading corrections into the leading ones. I checked 
that the trick does really help to improve the agreement in the first order between 4°' 7 ^ and 
the full 6\. But, as far as the TWIST experiment is interested in the region x > 0.3 (the 
event distribution is peaked at large x- values in any case), I don't apply the modification 
here. 

Let us define the relative contributions of pair corrections as 

5 (o,ee-) = V(a_\ 2 2/f NS) (x) +3/f S \x) ±cP,(2 g r S \x) + 3g?> 9 \x)) 



Q \2tT J /Born (a?) ± cP^Bc 

(l,ee) _ T ( a \ 2 ^ B ( x ) ± cP n92 ( x ) 



27T J ,- = n^s int /Born(x) ± cP^g Bo m( 



X j 



For two values of c, they are given in Table ^. One can see that the next -to-leading pair 
corrections have the same order of magnitude as the leading ones. This feature has been 
observed earlier in pair corrections to other processes |35|]. The functions, which describe 
the LL pair effect, have numerically small coefficients and are less divergent (at x — ► 1 and 
x — > 0) than the NLL ones. This means that the calculations of the non-logarithmic terms 
in the second order pair corrections is desirable, although the pair corrections are typically 
less than the photonic ones (at the same order in a). 

The combined effect of the discussed above higher order corrections to fj,~ decay spec- 
trum is shown in Fig. 

5 h . G . = 4°- 7) + 4 1 ") + + 4°' ee_) + + 4°' eg) + ££ p . (5.26) 

One can see from the Figure that the typical effect is of the order of 5 • 10 -4 for the 
intermediate range of x, and the effect becomes even larger at the two ends of the energy 
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c = 1 


c = -1 


X 


-i r\A A0.ee) 

10 4 <5 2 


lO 4 ^ ; 


-i r\d A0.ee) 

10 4 5g ' ' 


1 r\A A0.ee) 

10 4 <5 2 ' ; 


-i Al.ee) 

10 4 <5 2 ' ; 


-i r\d e(0,ee) 

10 4 5g ' ' 


0.05 


548.14 


-661.94 


-2.19 


2132.14 


-2848.25 


-9.12 


0.1 


66.02 


-68.36 


-0.18 


241.62 


-274.34 


-0.76 


0.2 


9.55 


—9.99 


n nr 

—0.05 


29.25 


on /^o 

—29.68 


—0.04 


0.3 


3.43 


-3.58 


-0.06 


9.59 


-9.77 


-0.05 


0.5 


0.49 


-0.71 


-0.06 


2.54 


-2.99 


-0.08 


0.7 


-0.86 


0.76 


-0.02 


0.46 


-1.03 


-0.08 


0.9 


-2.44 


2.95 


0.05 


-1.36 


0.99 


0.00 


0.99 


-4.95 


8.24 


-0.73 


-3.92 


5.47 


0.33 


0.999 


-7.32 


15.41 


-18.18 


-6.29 


11.76 


0.92 



Table 2: Pair corrections to \i decay spectrum versus x for c — ±1 and = 1. 
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Figure 1: The relative effect of higher order corrections versus electron energy fraction for different 
angles. 

spectrum. So the corrections under consideration are really important for the modern and 
future experimental studies of the muon decay spectrum, where elaborated technique and 
high statistics allow to reduce the experimental errors to the level of 1 • 10 -4 and better. 

6. Conclusions 

To estimate the theoretical uncertainty in the description of the polarized muon decay 
spectrum by Eqs. ( |2.1| ) and ( |5.1| ), we should consider the contributions, which have been 
omitted in the present calculations. They are: the O (a 2 ) order terms, which are not 
enhanced by the large logarithm; sub-leading contributions in the third order O (a 3 L m ) , 
m < 2; and all the leading and sub-leading effects in the forth and higher orders {O [a n L m ) 
where n > 4, m < n) except those ones, which are taken into account by exponentiation. 
The possible contribution to the uncertainty from strong interactions is negligible in our 
case, since it is suppressed at least by (a/V) 2 , and the lowest order the contribution of 
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hadronic virtual pairs was found in Ref. to be small itself. 

An estimate of the omitted contributions by a simple counting of powers of the fine- 
structure constant and the large logarithm is not very safe, because there could be some 
extra enhancement factors, like numerically large constant coefficients or powers of ln(l— x) 
(the latter is partially taken into account by exponentiation). I suggest to estimate the 
omitted terms by a linear extrapolation of the known expansions in a and L. Namely, 

xi 1 ) jc(°) 

x(2) _ r(l) 6 2 r(l) _ r(l) <>3 AO) _ .(0) $3 (R X 

d 2 T(oJ> ^3 ~<>2 -(oj. d4 d3 7(0)" ' ^ 

o 2 o 2 o 2 

(2) 

where denotes the contribution of the second order terms, which are not enhanced 
by any large logarithm; 5% is the third order next-to-leading correction (which can be 
calculated using the fragmentation function method described above); 5^ stands for the 
fourth order LL effect. This estimate of the theoretical error, can be applied to any par- 
ticular set of experimental conditions to derive the actual uncertainty. In principle, the 
latter can depend on various cuts and details of particle registration and event selection 
(see discussion in Ref. @l). The approach (|6.1[) to estimate the uncertainty in the muon 
decay spectrum description works well for the main part of the kinematical domain. But, 
if one studies separately the extreme region x <C 1 (or l-x< 1), a special investigation 
of the convergence properties of our perturbative expansions in a and L should be per- 
formed. Evaluation of the uncertainties for any concrete experiment can be done using the 
analytical results and applying specific conditions of particle registration and data fitting. 
The new results of Ref. [18] and the present paper reduce the theoretical uncertainty in 



the description of the polarized muon decay spectrum. For the quasi-realistic experimental 
setup described in Ref. |Q, we can obtain now about 1.5 times better precision, so that, for 
instance, the theoretical uncertainty for the Michel parameter p becomes 2 • 10 -4 instead of 
3-10 -4 obtained in Ref. [fL4f |. Nevertheless, this is still worse than the experimental precision 
1 • 10 -4 planned at TWZST ||, Assuming that a theoretical precision of about one 
third (or less) of the experimental one would not spoil results of an experiment, we see a 
challenge for further investigations. First of all, a calculation of the O (a 2 ) contributions, 
which are not enhanced by the large logarithm, is required to ameliorate the theoretical 
precision. This calculation is difficult, but possible by means of the standard methods. 

The formulae for higher order corrections (with simple substitutions) are valid also for 
the decays of r-lepton: r — ► \iv T v^ and r — > ev T v e . 
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